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We analyze the harvesting of entanglement and classical correlations from the quantum vacuum 
to particle detectors. We assess the impact on the detectors’ harvesting ability of the spacetime 
dimensionality, the suddenness of the detectors’ switching, their physical size and their internal 
energy structure. Our study reveals several interesting dependences on these parameters that can 
be used to optimize the harvesting of classical and quantum correlations. Furthermore, we find 
that, contrary to previous belief, smooth switching is much more efficient than sudden switching in 
order to harvest vacuum entanglement, especially when the detectors remain spacelike separated. 
Additionally, we show that the reported phenomenology of spacelike entanglement harvesting is not 
altered by subleading-order perturbative corrections. 


I. INTRODUCTION 

It has been known for a long time that the vacuum 
state of a free quantum field contains correlations be¬ 
tween time- and spacelike-separated regions [U [2] . Be¬ 
sides the remarkable fundamental interest from the point 
of view of quantum foundations, the existence of this vac¬ 
uum entanglement is a key ingredient in very interesting, 
recently-discovered phenomena such as Masahiro Hotta’s 
quantum energy teleportation Eg. It is also at the core 
of long open problems such as the black hole information 
loss problem [5] and some of its proposed tentative so¬ 
lutions, such as the so-called “black hole firewalls”, and 
black hole complementarity [MSl. 

A perhaps more surprising result is that this vac¬ 
uum entanglement can be extracted from the field to 
Unruh-DeWitt particle detectors [9] that couple to the 
field locally even when the two detectors are spacelike- 
separated, as pointed out, first by Valentini m, and later 
by Reznik mma. This phenomenon has become known 
as entanglement harvesting m- 

Since Unruh-DeWitt detectors can, in some regimes, 
be a good approximation to the light-matter interaction 
[HlIIS], these pioneering results may imply that it is pos¬ 
sible to extract entanglement from the electromagnetic 
vacuum to atomic qubits, where it could be used as a 
resource, although, for this purpose, one has to be care¬ 
ful with the impact of time synchronization on entangle¬ 
ment harvesting [16]. Indeed, it has been proved that it 
is possible to devise quantum optical setups where entan¬ 
glement can be sustainably and reliably extracted from 
a quantum field and distilled into Bell pairs, that can be 
later used as a resource for quantum information tasks. 
This technique is known as entanglement farming m- 
Moreover, there have been several exploratory works on 
the experimental feasibility of timelike and spacelike en¬ 
tanglement harvesting in atomic physics and supercon¬ 
ducting circuits mmo]. 

Entanglement harvesting was proven, by Ver Steeg and 
Menicucci, to be sensitive to the structure of the back¬ 
ground spacetime in which it is performed [21] . In partic¬ 


ular, they proved that entanglement harvesting can dis¬ 
tinguish between a thermal background and the Gibbous- 
Hawking radiation background of an expanding universe 
[2IH23], and it is also sensitive to the topology of space- 
time [23] . 

Not only that, entanglement harvesting has been 
proven to be very sensitive to the state of motion of the 
detectors, and the boundary conditions on the field on 
which is performed. This has led to proposals of appli¬ 
cations in metrology such as range-finding m and even 
quantum seismology [25] . 

As a fundamental phenomenon, entanglement harvest¬ 
ing is therefore relatively well understood. However, little 
is still known about how it is affected by (and possibly op¬ 
timized over) variations of the specific parameters of the 
setup, such as how fast the detectors are switched on, the 
dimension of spacetime, the physical size of the detectors 
or the nature of their internal degrees of freedom. This is 
particularly important in the case of spacelike entangle¬ 
ment harvesting, which would constitute a direct proof of 
the existence of vacuum correlations. We can find several 
hypotheses and intuitions about some of these aspects 
in the literature. Eor example, in one of the original 
papers by Reznik et al [12], in which they maximize 
entanglement harvesting using very fast-varying super- 
oscillatory switching functions, or in [26] (in the context 
of harmonic oscillator-based non-perturbative methods 
for particle detectors Eaui]), where it was speculated 
that a sudden switching might be more efficient than a 
smooth one to harvest entanglement. 

In this paper we present a thorough study of both en¬ 
tanglement harvesting and the harvesting of classical cor¬ 
relations, and how they are affected by the dimension of 
spacetime, the physical size of the detectors, their inter¬ 
nal energy structure and the smoothness of the switch¬ 
ing of their coupling to the field. Remarkably, and con¬ 
trary to previous belief, we find that smoother switchings 
are much more efficient than sudden switchings in order 
to harvest vacuum entanglement. Namely, we find that 
while for a smooth Gaussian switching it is always pos¬ 
sible to choose detector setups that allow for spacelike 
entanglement harvesting, this is not the case for sudden 
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switchings. We trace back this result to the fact that 
sudden switchings increase the amount of local noise that 
the particle detectors experience [28] , which hinders their 
ability to harvest vacuum entanglement. 

We show that entanglement harvesting is rather in¬ 
sensitive to the dimensionality of spacetime, but this is 
not the case for the harvesting of classical correlations. 
Namely, for a 3+1-dimensional spacetime, mutual infor¬ 
mation is more efficiently harvested from the vacuum 
than in 1+1 dimensions when the detectors are in light- 
like contact. On the other hand, reducing the dimension¬ 
ality of spacetime improves the ability of the detectors to 
harvest correlations when they are spacelike separated. 

We also show that finite-size, but small detectors (as 
compared to their interaction time with the field) do not 
behave in a fundamentally different way to point like de¬ 
tectors, in the regimes where the point like approximation 
is not ill defined. When the size of the detectors is in¬ 
creased and it becomes comparable to their interaction 
times, larger detectors are much less efficient to harvest 
entanglement than smaller ones. 

As for the dependence on the detectors’ energy gap, 
we show that the situation is radically different in the 
cases of sudden and smooth switching. For the latter 
it is always possible to tune the detectors’ energy gap 
in order to harvest spacelike entanglement for a given 
setup, whereas for a very fast switching it is generally 
not possible to do so. 

Finally, we have also analyzed vacuum entanglement 
harvesting at higher orders in perturbation theory, show¬ 
ing that going beyond leading order does not reveal new 
phenomenology. Therefore a leading-order perturbative 
approximation is generally enough to identify the regimes 
in which Unruh-DeWitt detectors can harvest quantum 
entanglement from the field vacuum. 


II. SETUP 

We will model two particle detectors (A and B) with 
the well-known Unruh-DeWitt model [9]. Although sim¬ 
ple, this detector model comprises most of the fundamen¬ 
tal features of the light-matter interaction when there is 
no exchange of angular momentum [HIITS]. The Unruh- 
DeWitt detectors (from now on referred to as the ‘atoms’ 
or ‘detectors’) interact with a background scalar field via 
the following Hamiltonian in the interaction picture, 

^ f F^{x - (1) 

r/e{A,B} 

where \y is the overall coupling strength, is the 

monopole moment of each detector 

fj'uit) = (T+ (2) 

(cr^ are SU(2) ladder operators), Fy{x) are the spatial 
smearing functions of each detector and Xj, their respec¬ 
tive center-of-mass positions. Xzy(t) is detector z/’s switch¬ 
ing function, which controls the interaction time and the 


coupling strength of each atom with the field. For our 
purposes we are going to consider switching functions 
that are strongly suppressed outside of finite time inter¬ 
vals (so as to have finite-duration interactions). 

Typically, if we think of the UDW model as a model 
of the light-matter interaction, the spatial support of 
the atom can be associated with the spatial probabil¬ 
ity profile of the atomic wavefunctions m- In an n+1- 
dimensional fiat spacetime, the scalar field can be ex¬ 
panded in terms of plane-wave modes in the following 
way 


(f>{x,t) = 


/ 


d"fe 


v/(27r)"2|fc| L 


+ H.c. 


( 3 ) 


The creation and annihilation operators ap. and satisfy 
canonical commutation relations [afe,at] ~ — k'). 

The integral with respect to x can be easily performed 
yielding the Fourier transform of the spatial profile. 


Hi = f 

V 


d"fc 

\/2|fc| 


( 4 ) 


where the form factor of the atoms is defined as 


AW=y=/d”-X(x)e'‘-. (5) 


From now on, we will consider that both detectors will 
have the same shape and their spatial profile will be a 
real function, therefore Fj,{k) = Fj,{—k) = F{k). 

The time evolution generated by § can be obtained 
perturbatively through a Dyson expansion of the time 
evolution operator: 


/ CxD nOO nt 

dt dt'H{t)H{t') + ... (6) 

-oo J—oo J— oo 


t/(l) 


Z7(2) 


If the initial state of the detectors-field system is po, 
the evolved state will be given by p = U. Let us 
define the notation O (X^) to represent terms that are 
proportional to V /c, / such that k X I = n. We will 

denote the O (A^+-^ ) as a sum of terms of the form pertur¬ 
bative contributions to the time-evolved density matrix 
as 


(7) 


Therefore we can write the time-evolved density matrix 
as a sum of terms of the form of 0 : 


p = + +p(0.1) +^(2.0) +^(0,2) ^^(1,1) +^(1.2) + _ . _ 


(8) 

Since we are going to analyze entanglement and cor¬ 
relations harvesting from the vacuum, we consider that 
the initial state of the detectors-field system is 


Po = |0)(0| ® PABfi, 


( 9 ) 
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where |0) is the vacuum state of the scalar field and Pab,q 
is the initial state of the detectors. We will be interested 
in the partial state of the detectors after their interaction 
with the field, which is given by 


PAB = TV^([//9oC/t). (10) 


This means that the non-diagonal terms in the field pro¬ 
duced by time evolution will be of no relevance for our 
purposes. In particular, any contribution for which 
the parities of i and j are different (i.e, for all the con¬ 
tributions with odd powers of the coupling strength Xj,) 
will give a zero contribution to the detectors’ final state 
(10), as long as the initial state of the field is diagonal in 
the Fock basis (as it is the case of the vacuum). 


Consequently, in the perturbative expansion of Pab^ 
the first order correction is trivially zero for the rea¬ 
son explained above. To leading order in the coupling 
strength, the two detectors’ time-evolved density matrix 
is given by 


pAB — pAB,0 + P^aB + P^AB + P^AB + (H) 

where = Tr^ [p'^'d) j, 

We are going to consider the case in which both detec¬ 
tors are in their respective ground states 


pAB,o = \gA){gA \ C) \gB){gB\- (12) 


From (11), pab takes the following matrix representa¬ 
tion 


PAB 


/I — Caa — ^BB 0 0 Al*\ 

0 Caa P^ab 0 I 

0 Cba Cbb 0 I 

V M 0 0 0 / 


in the basis 


^0{Xt) 

(13) 


{|Pa) C) Ipb), \^ a ) C) Ipb), Ipa) C) le^), |e^) (g) le^)} • 

(14) 

The explicit expressions of and M can be easily 
obtained from Eq. after substituting 0, 0 and 0, 
yielding 


= j d^kL^{k)B{k)*, (15) 

M=[d^kM{k), (16) 


where L^{k) and M{k) are 


L —Ap 




F{k) 




/ oo 

-OO 


(17) 


M{k) = - 


[HkW 


2\k\ 

/ OO nt\ 

dh / 

-OO J —oo 

+XB(ii)XA(i2)e'(“®*^+“-^‘=) 


(18) 


The expressions above are rather general and can be 
easily particularized to any switching and spatial profiles 
for any dimension. 

We would like to analyze under which set of general 
conditions it is possible to harvest classical correlations 
and entanglement from the field to the detectors. With 
this aim, we will consider different switching modalities 
(sudden versus Gaussian), different characteristic detec¬ 
tor sizes (pointlike versus non-negligible Gaussian smear¬ 
ing), different spacetime dimensions (1+1-dimensional 
—as in long wave guides or optical fibers- versus 3+1- 
dimensional —as in free space-) and a range of different 
detector internal energy scales. 

The first step is to evaluate the integrals ( pT[ ) and (18) 
for the different cases that we will consider. Let us first 
focus on the spatial profile. We will choose the following 
Gaussian smearing 


F{x) = 


1 




or.2 /^2 

^-X l(J 


(19) 


which in turn enters equations 0 and ( p!^ via its 
Fourier transform 


F(k) = 






( 20 ) 


It will be relevant to consider the limit where the detector 
is pointlike localized in space (a ^ 0). Namely 


F{x) = ^ F{k) = 


1 




( 21 ) 


A. 3+1 dimensions 

We will first consider the case of three spatial dimen¬ 
sions. We will study the following scenarios 

Al. Gaussian switching and Gaussian spatial smearing. 

A2. Gaussian switching and pointlike detectors. 

A3. Sudden switching and Gaussian spatial smearing. 

A4. Sudden switching and (almost) pointlike detectors. 
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Al. Gaussian switching functions and Gaussian smearing 

Let us first consider the case in which the detectors 
are turned on in a smooth manner, following a Gaussian 
profile 


Xv{t) 




( 22 ) 


With this switching function, all time integrals in 
( [T7[ ) and ( p! 8 | ) admit analytic closed forms. We will 
assume that both detectors have the same energy gap 
Vt = VLa = and that both couple to the field with the 
same strength A = = A^. This allows us to scale all 

the parameters in the system relative to the characteris¬ 
tic timescale of the switching function T, as suggested in 
[niin]. Defining the dimensionless magnitudes a = DT, 
/3/x = S = (JjT^ K, = kT and = t^/T (a sum¬ 

mary of all the dimensionless parameters used through¬ 
out the paper can be found in Table |I|, ( [Tt] ) and ( p^ can 
be recast as 

Gi{K,Tf,), (23) 

2 |k| {2Trf 

_ 1 ^ 2(52 

M(fc) = - - oG2iK). (24) 

2\k,\ ( 27 r) 

where, taking = tijT as dimensionless integration vari¬ 
ables, we can write 


L^(k) =AT 1/2 


Gi{K,Tf,)=T (25) 

J —OO 

/ CO nTi 

-OO J —CO 

X , (26) 


where 7 = /T is the normalized separation be¬ 

tween the switching functions’ centers. 

The integral ( [^ can be readily analytically evaluated. 
Obtaining a closed form for (26) is more involved, but it 
can be accomplished via parametric differentiation un¬ 
der the integral sign and solving the resulting differential 
equation (see details in Appendix [A|). The result in both 
cases is 


Dimensionless 

variable 

Expression 

Physical meaning 

Oi 

QT 

Energy gap 

/3m 

Xf,lT 

Detectors’ positions 

P 

d/T 

Spatial distance 

7 

A/T 

Time delay 

(5 

o/T 

Detectors’ size 

K/, T] 

kT, qT 

Momenta 

T 

t/T 

Time parameter 


TABLE 1. Collection of all the dimensionless quantities that 
are used throughout this paper. Notice that d — \xb — ^ca|, 
A = tB-tA and (3 = \f3B - /3a\. 


(16), yielding (see Appendix [A|) 


^AA = 


47r2T2 




pCO 

- / d|K| |k| e~^'^''^^Gl{K,TA)Gl{K,TA) 

Jo 


87r(l + (52) 


VTVP 


(29) 


Gab = 


A 2 


pCO 

J d|K| sin(|K;|/?)e“5''^'’^Gi(K,0)Gi(K,7) 


47r2T2/3, 

' 87^/3 \/rT^ 






\M\ = 


A 2 


47r2r2/3 


V 2 VTTP 

nCO 

/ d|K| sm{/3\K\) e~G 2 {k) 
Jo 


(30) 


A^e 2 “ 


X [E{k,A + E{k,-A] , 


pCO 

/ d|K| sin(/?|K|)e-ni+^')«' 
Jo 


(31) 


Gi{k,t^) = (27) 

G2{k) = ^ _ 

where, for simplicity, we define 




1 — erf 


7 + i|/«|A 

^/2 J 


(28) 


Using (23), (24) and (27) we can evaluate (15) and 


where we have defined p = 1/3^ — /3a | = d/T where 
d = \xb — xa \ is the distance between the detectors’ cen¬ 
ters of mass, and erfc(x) = 1 — erf(x) is the complemen¬ 
tary error function. Note that, with the assumptions 
made, Cbb = Caa- 

This expression is general for any separation between 
the centers of the Gaussians 7 = (tB—'^j^jT. We 
note that the integral over |k| admits an approximate 
analytic closed form when 7 is large enough to neglect 
the overlap between the two Gaussian switchings. To 
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see this, we start from (18) and notice that Xu{t) are are in perfect overlap (i.e, 7 = 0). In this case 


Gaussian functions with standard deviation st = Tj\f2. 
When the centers of the two Gaussians are separated by 
A = tB-tA> 7T/^/2 (i.e. more than seven times the 
standard deviation), the two Gaussians effectively do not 
overlap. Note that, for A > 7T/v^, the overlap between 
the Gaussian switchings is suppressed by a factor smaller 
than ^ 10“^^, this is, any effect of that overlap 

will be suppressed at least 10“^^ times as compared to 
the contributions of the vicinity of the Gaussian max¬ 
ima, rendering that overlap region completely negligible 
for our purposes. Assuming without loss of generality 
that the detector A is switched on before the detector B, 
in this case the first contribution of (18) is very approx¬ 


imately zero, and \JVl\ ^ |Alnon|, where 


iMnonI = 




8v^/3\/r+72 

(/3-7)^ r 
e 2(1+52) 


X 


(^+ 7 )^ 

_e 2(1+52) 


1 -1- erf 
1 — erf 1 1 


A P-i A 

V\/2+rT+/ 


/3 + 7 


^/2^/lTP J. 


(32) 


A2. Gaussian switching functions and pointlike detectors 
For pointlike detectors we take the limit (5 ^ 0 in Eqs. 


(29), (30), (31) and (32), yielding 


C 


AA 


}?e 2^ 
SttP 


2 — 02^ aerie 


, (33) 


^AB = 




sV^/s 

(/3-7-ia)2 


X < e 


erfc i 


. /3 — 7 — ia 


V2 


(7 + /3 + ia)2 / fd J^ia 

—e 2 erfc —1 


\M\ = 


X^e 2^ 


V2 

poo 

/ dj^l sin(;d|K|) e~^ 

Jo 


(34) 


SttP 

X [-^(k, 7 ) +-^(k,- 7)] 


(35) 


IMnonl = 


A^e 2“^ 

8+^/3 

(/3-7)^ 


1 + erf 




V V2 


(P±jy 


—e 


1 — erf i 


■P+i 


V2 


(36) 


It is worth mentioning that in this case there is another 
situation in which A4 has an analytic closed form, which 
is when the two identical detectors’ switching functions 


Me 


y2g-|(a2+/32) 


erfi ( ^ 1 - i 


(37) 


where erfi(x) = —ierf(ix) is the imaginary error function. 

AS. Sudden switching functions and Gaussian smearing 

Let us now consider that the detector u is switched on 
in an abrupt manner at and switched off in the same 
way at Namely, the switching functions in 0 and 
(18) are now 


Xuit) = 


1 if T°" < t < Tf, 
0 otherwise. 


(38) 


For this switching function, the time integrals also ad¬ 
mit closed-form expressions. From (pTl and (18) we get 


L^(fe) =xB/^ 




2\k\ {271) 




(39) 


e 2^ 


rS2{K), (40) 


M{k) =_ 

2 |k|( 2+" 

where, similar to the Gaussian case, T = r°® — T°" is 
the timescale of the interaction, and 5 'i(k, r^) and S2{k) 
are defined by 

5i(«,t^) = T / " (41) 

J^on 

/ OO pti 

-OO J —OO 


X [XA{ti)XBit2) + XBitl)XAit2)] 


(42) 


As in the Gaussian case, (4n can be computed in a 


straightforward manner, while the case of (42) needs a 
careful analysis (see Appendix]^. These integrals yield 

iT 

S.(-=,rJ=--, (43) 

and S 2 {i^) simplifies to 

5^2 (^) = ^2non(^) (44) 

^ rr.2 _LMT_ ll 

-L O n n 

— a^ 

when Tg" > T/ (no overlap between the switching func¬ 
tions), and to 

(45) 


^2(^) ^2over(^) 

rj-i2 


+e 


gi(7+l)(«-|«^l) ^g2i7|/«| _ gi(a+|/«|)^ 

(e2i«7 _ g2ia^ 


,i(7+l)(a-|«^|) _ i7(a-k|) _ rjo (^p2ia7 _ 2ia\ 
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when the detectors’ switching functions overlap 
< Tf. 

Inserting 


and (16) take the form 


(44) and (45) into (39) and (40), (15) 


same as in the 3+1-dimensional scenario that we just 
computed. Therefore, in the 1+1-dimensional case, 
and \Ai\ take the form 


\ 2 roo 

Caa=^ 

TT^ 


/ dM® 

Jo ( K 




(l«l +«)" 



g Ve-s"" 

‘ 1 

Caa = . 

2(“+ l^^l) 

4 

5 

nOO 


A 1 -i 

/ dK -—re 2 

J-oo W 




(46) 


+ 


Cab 


47r^/3 


e“ 

Jo ( K 




Ja i^i 


-Ap-h'^ (i+<5 


(50) 


(|k| +a)' 


-sin(/3|K|) 


Cab = 




Inon I — 




f . e 2^ ^ 

j„ 


(47) 

X 


_ 


r-A IK.P 

/ dK — e^ 

J-oo 1^1 

J 


47r^/3 

X - l) - l) 


pi/t/3 


(48) 


Inon I — 




_gi + l7g-|/^ (1+^ )g- + + 


+ 1 
^-A 


, ( 51 ) 


\.hv ^ 

b\ 


\M, 


A^ 


47r^/3 


POO 

/ d|K| sin(^|K;|)e“5'5^«''5'2over(K) 


Jo 

Ja I^I 


(49) 


where Alover and Alnon correspond to the values of A1 
when the detectors’ switching functions overlap and when 
they do not, respectively. 


A4- Sudden switching and quasi-pointlike detectors 

It is well known that a sudden switching for a point¬ 
like detector in 3+1 dimensions leads to ultraviolet diver¬ 
gences in the response of a particle detector [28]. How¬ 
ever we can consider detectors whose spatial smearing is 
much smaller than the duration of the interaction (i.e. 
cr/T 1) as an effective pointlike detector that can be 
switched abruptly in a 3+1 dimensional scenario. 


B. 1+1 dimensions 

As in the 3+1-dimensional case, we will explore the 
following different switching and spatial profile configu¬ 
rations for the detectors: 

Bl. Gaussian switching and Gaussian spatial smearing. 
B2. Gaussian switching and pointlike detectors. 

B3. Sudden switching and Gaussian spatial smearing. 
B4. Sudden switching and pointlike detectors. 

Bl. Gaussian switching functions and Gaussian smearing 


\M, 


X^e 


L 




■"/a 


— oo 

oo gi/^/3 

A l« 


K\ 


(52) 


£;(k,7)+£'(k,-7) 


+/ dK—e-^'* (i+'5 ) 


-B(k,7) + -E(«,- 7) 

(53) 


where the notation Alover and Alnon is understood in 
the same way as in section IIA, and A is an infrared 


cutoff which regularizes the well-known logarithmic IR 
divergences of the 1+1-dimensional case. This kind of 
IR regularization is common in the literature and can be 
thought, for instance, as the length scale of a very long 
optical cavity or a periodic optical fiber, or the charac¬ 
teristic radius of a cylinder spacetime topology. 

Notice that, as opposed to the 3+1-dimensional case 
(where the field has units of [(j)] = T~^ and A is dimen¬ 
sionless), in 1+1 dimensions the field is unit less, hence 


A has units of T. Therefore, in Eqs. (50), (51), (52) 


and ( [53| ) we have already used a dimension! 
strength defined as A = AT. 


ess coup 


ing 


B2. Gaussian switching functions and pointlike detectors 

As in the 3+1-dimensional scenario, the pointlike case 
corresponds to the limit (5 —> 0 in Eqs. (50), (51), (52) 


and (53), yielding 


When we particularize (15) and (16) to 1+1 dimen¬ 
sions, the time integrals in (17) and (18) are exactly the 


Caa=^ 


L 


-Ag-|(a+ + |)2 roo ^_i(a+ + |) 


dK, 

— 00 


■/ 

Ja 


dhc- 


, (54) 
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= 


A2e-5«"ei«7 


|*^non| — 


|Xc 


4 


POO 

Ak/3 

+ / dfC-r-r 

Ja 

kl 

Ve-A" 

r 

4 

L 

POO 

AkI: 

+ / d/t —- 

Ja 

\K\ 

Ve-i 


8 


/ 


—A 11^13 

k 


yg-5«^g-|«l« 


dK—e‘kl7e-i«" 

k 


/d/t——e ^ 

J-00 kl 

E{k,')) + E{k,-P 

+/dK——e 

Ja kl 

£;(k,7) + £;(k,-7) 


(55) 


(56) 


(57) 


|•^^over| —“j 

47r 


/ d/D^e 2 S2over{l^) 

J—oo \^\ 


roo 

+ / dK-p-pe “2 S2over{K) 

Ja 




(61) 


B 4 . Sudden switching and pointlike detectors 

In 1+1 dimensions the point like limit of the previous 
case is not divergent, and we can therefore safely take 
the limit J ^ 0 in Eqs. (58), (59), (60) and (61), which 
results in 


-4 / 


dhv ■ 


1 


-00 |/^| (|/^| + a ) 


h« + kl) 


i 


dK ■ 


A \k\ + a) 


sm 


:(«+ kl) 


(62) 


B3. Sudden switching functions and Gaussian smearing 


In this case the time integrals are again (41) and (42). 
Inserting them into 0 and ( p+l ) along with the smearing 
function (20), (H!^ and (H!^ now read: 


^-4 / 


dhv ■ 




—00 

CXD 


dK - 


kl(kl+ + 

e-5«"'5" 


■ sm 


24+ l«^l) 


/A pc (pc +(a) 


-((7 + kl) 




/ 


-A IK3 

dK 


k\ (kl + of 

X g-i(7+l)(a+|K|) ^gi(a+|K|) _ 

+ / d/D-^ 

X g-i(7+l)(a+ + l) ^gi(a+ + |) _ 




too |/t| (|/v|^ — aA 

X ^ei(«-|«l) _ ^ei(«+kl) _ 
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We have now all the ingredients needed to study both 
entanglement and general correlations from a quantum 
scalar field in 3+1 and 1+1 flat spacetimes for the differ¬ 
ent switching functions and spatial smearings considered. 


III. ENTANGLEMENT HARVESTING 

We are interested in quantifying the entanglement ac¬ 
quired by particle detectors after their interaction with 
the vacuum state of the field. As we will see, factors 
such as the dimensionality of spacetime, the nature of the 
switching functions, the smearing of the detectors and 























































their internal energy structure have a dramatic impact 
on their ability to harvest field vacuum entanglement. 

We will use negativity [29] to quantify entanglement. 
For a two-qubit system, the negativity (defined as the 
sum of the negative eigenvalues of the partially trans¬ 
posed density matrix) is an entanglement monotone 
which only vanishes for separable states [30l |3T] . 


A. Entanglement harvesting to second order in 
perturbation theory 


To second order in perturbation theory there is only 
one eigenvalue of the partial transpose of (13) that can 
be negative 




— \J {^AA — ^Bb)^^^ \MY 


( 66 ) 

Note that a naive inspection of the partial transpose of 
would have produced the apparently always negative 


eigenvalue E 2 = — {jCab^- However note that {jCab^ is 
O (A^), and thus £’2 =0 + 0 (A^)- As we will discuss 
later, to find the correct form of £2 the whole fourth 
order correction to the density matrix must be computed, 
resulting in £2 not being negative in general. Therefore, 
we define the following negativity estimator: 

= - £1 


^AA + ^BB — \J {^AA — ^Bb)^ + 4 \M\^ 


(67) 


so that the negativity, to second order in the coupling 
strength, is equal to A7 = max(0, A7^^^). 

This form of the negativity estimator can be further 
simplified when it is assumed that the two detectors are 
identical and switched on for the same amount of time 
(but with a time delay between them), and therefore 
= Caa = ^bB‘ Under this assumption, at leading 
order in perturbation theory the previous formula can be 
written as 


=M-C^ 


( 68 ) 


This expression shows very intuitively that for the state 
to be entangled, the nonlocal term A4 has to be larger 
than the local terms [12]. We can now evaluate the 
impact of the smoothness of the switching on the ability 
of the detectors to harvest entanglement from the vac¬ 
uum. In view of (68), it is not straightforward what the 


impact of the switching may be: a sudden switching may, 
on the one hand, increase the local noise (the particle 
count of the individual detectors, Caa and Cbb)^ but on 
the other hand it is not a priori clear what would be the 
impact of the switching on the nonlocal terms A4. There¬ 
fore we will need to compare the two scenarios (smooth 
vs. sudden switching) to confirm if the extra noise in¬ 
troduced by a sudden switching hinders the detectors’ 
ability to harvest entanglement. 


In Fig. 12 we present some examples of the negativity 
in all the cases defined in sections III Al and IIIBI For the 
binary plots showing the parameter region where entan¬ 
glement harvesting is possible (a.l, b.l, c.l, d.l, e.l, f.l, 
g.l, h.l), we have chosen to fix A/T = 0 (no delay be¬ 
tween the switchings) for the sudden switching cases be¬ 
cause this maximizes the parameter region where entan¬ 
glement harvesting is possible. This is because abruptly 
switched detectors quickly lose their ability to extract 
vacuum entanglement as A/T grows, even in lightlike 
connection. For the Gaussian switching this is not the 
case, and therefore we have chosen to show the param¬ 
eter regions of entanglement harvesting for A/T = 3 in 
order to better display the features of the parameter de¬ 
pendence. 

The plots (a.2, b.2, c.2, d.2, e.2, f.2, g.2, h.2) show 
the variation in the negativity estimator with 

the distance between the detectors (d/T) and their 
switching delay (A/T) for a representative value of 
OT = 7. Notice that, since the negativity is defined as 
A7 = max(A7^^\0) + O(A^), in the cases where the plots 
display negative values there is no entanglement harvest¬ 
ing at all. It is still illustrative in those cases (all of them 
sudden switching) to show how far from zero the nega¬ 
tivity estimator would be. 

First, we note that the entanglement peaks in the re¬ 
gion where there is lightlike contact between the detec¬ 
tors (in between of the two white, dashed lines). In this 
situation both detectors can exchange real quanta via 
the background field, resulting in a peak in the nega¬ 
tivity that will be also visible in the study of the to¬ 
tal correlations. In the majority of the cases negativity 
peaks when there is full lightlike contact (exactly at the 
exact center of the region inside the two white lines in 
Fig. [^. However, notice that there are two different 
trends: the entanglement increases when the detectors 
are light like connected, but decreases in a nontrivial way 
with the spatiotemporal distance, thus it is not surprising 
that in some cases (for instance, for the I+I-dimensional 
cases with sudden switching functions) the maxima of the 
negativity estimator do not appear at the center of the 
light like region, although these maxima are still always 
somewhere in the region of light contact. 

Another clear feature that arises from our analysis 
is that, perhaps contrary to intuition, the detectors 
switched on in a smooth, more adiabatic manner are ca¬ 
pable of harvesting entanglement in a much more efficient 
way than suddenly switched detectors, both when the de¬ 
tectors are in causal contact and when they are spacelike 
separated. 

Let us focus now on the interesting case of spacelike 
entanglement harvesting. Since the detectors cannot ex¬ 
change information, spacelike entanglement is the result 
of the harvesting of only preexisting correlations in the 
quantum vacuum. 

We see from Fig. [^that for smooth Gaussian switch¬ 
ing it is always possible to obtain spacelike entanglement 
harvesting if the internal energy gap of the detectors fl is 
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FIG. 1. Negativity for all the cases described in sections IIA and in the following order: (a) three-dimensional (3D) 
Gaussian switching and Gaussian spatial profile; (b) 3D Gaussian switching and pointlike detectors; (c) 3D sudden switching 
and Gaussian spatial profile; (d) 3D sudden switching and near-pointlike detectors (cr/T ^ 1); (e) one-dimensional (ID) 
Gaussian switching and Gaussian spatial profile; (f) ID Gaussian switching and pointlike detectors; (g) ID sudden switching 
and Gaussian spatial profile; (h) ID sudden switching and pointlike detectors. The plots a.l), b.l), ..., h.l) (first and third 
columns) show in (dark) red the values of d/T and QT for which entanglement harvesting is possible, and in (light) grey the 
values for which there is no entanglement harvesting. The plots denoted by a.2), b.2), ..., h.2) (second and fourth columns) 
show the specific value of the negativity estimator as a function of the spatial separation of the detectors and the time 

delay between their switching functions for a given value of QT (recall that there is entanglement harvesting when > 0). 
In all plots the dashed lines represent the boundaries of the lightcone. In the second and fourth columns, the points to the right 
of the rightmost dark line represent spacelike separation. These boundaries are placed at A = d ± T in the cases of pointlike 
detectors which are switched suddenly, and at A = d± 7T/in the rest, which is a reasonable estimation given the discussion 
in section II, subsection Al. All cases of Gaussian spatial profile have cr/T = 1 except for the near-point like detectors (plots 
d.l and d.2), for which a/T = 0.01. We see that Gaussian switching always allows for spacelike entanglement harvesting for 
sufficiently large values of QT, while this is not the case for sudden switching. For the 1+1 dimensional scenarios the IR cutoff 
was taken AT = 0.001, much smaller than all the relevant frequency scales in the setup. 
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increased (consistent with the analysis in [l3]), but pay¬ 
ing the price that the larger the value of the smaller 
(but finite) the amount of entanglement harvested. In 
other words, for Gaussianly switched detectors, when 
the energy gap increases the total amount of entangle¬ 
ment decreases rapidly, but in exchange the region of 
nonzero negativity increases, eventually “leaking” arbi¬ 
trarily much into the spacelike-separation region. 

Strikingly, this is not the case for sudden switching, 
where the behaviour of negativity with the internal gap 
of the detectors is dramatically different to the Gaussian 
switching case. 

In particular, for sudden switching, negativity is no 
longer monotonic with the detectors’ energy gap and, 
what is more, the values of ft that allow for spacelike 
entanglement harvesting are severely limited. We found 
that there are even cases in which it is impossible to 
harvest entanglement with spacelike-separated suddenly- 
switched detectors for any value of Q within the range 
explored. Increasing values of ft do not seem to improve 
this situation (see Figs. [Il> Et. Ef); it is only possible 
for point like detectors in 1-i-l dimensions to find space¬ 
like entanglement harvesting with a sudden switching, 
and increasing ft actually reduces the spacelike harvest¬ 
ing ability even in this case (see Fig. Ef)- 

The conclusion we extract is that the increase of the 
local noise as a result of a sudden switching hinders our 
ability to harvest quantum entanglement from the field, 
regardless of any possible beneficial effect that the sud¬ 
den switching may have in the nonlocal A4, whether the 
detectors are or not in causal contact or the dimension 
of spacetime. 

Finally, let us consider how our ability to harvest en¬ 
tanglement from the vacuum varies with the size of the 
detectors. This behaviour is shown in Fig. First note 
how the behaviour of the detector quickly goes to the 
pointlike limit when a/T 1. 

For larger a we see that as the sizes of the detectors 
grow, they become less efficient to harvest vacuum en¬ 
tanglement. This is true as long as a < 6d, this is, when 
the two detector Gaussian smearings are further apart 
than about 8 standard deviations. Otherwise, if we were 
to allow the detectors to have a spatial overlap (a rather 
unphysical situation, however), this overlap enhances the 
ability of the detectors to get entangled. 


B. Entanglement harvesting beyond second order 
in perturbation theory 

Recall that in the analysis of the negativity that we 
carried out in section IIII Al we found that there was an 
eigenvalue of the partially transposed density matrix, E 2 , 
which was O(A^). This eigenvalue ‘naively’ appeared to 
be always negative when taking into account only the 
contributions coming from the square of second-order 
perturbative terms. However, to find the right fourth- 
order perturbative corrections to the negativity (and in 


particular, the right value of E 2 ) it is not enough to cut 
the perturbative expansion at second order as in (13). If 
we expand pab to O(A^) we obtain 


PAB 


/ Pll 

0 

0 

0 

7-2^2 

P 23 

0 

P 23 * 

P 33 


0 

0 


PlE 

0 

0 

Paa/ 


+ 0{\l), (69) 


where 


(Si + S 2 + S 3 ), 


(70) 


Pli =1 — {^AA + ^Bb) 

P^22 =^AA + ^ 1 , 

P^33 =^BB + ^2, 

P44 ^3 5 

pZ =M^ + T, 

P 23 =^AB + n. 


Here, all the Greek letters denote terms which are pro¬ 
portional to the fourth power of the coupling strength. 
The partial transpose of (69) has two potentially negative 
eigenvalues: 


El =- 

2 


(f^AA 


^BB — 


\/{^AA — f^BB)^ + 4 


+ Si + 


4Re(TA/l) 


4 \Mf‘ + {Caa — ^bb) 
{^AA — ^bb) (Si — S2) 


^ 4 \Mf + {Caa — ^bb)‘^ 
E2={Es-\Cab\^)eO{XI). 




(71) 

(72) 


El is just the eigenvalue found in the study of 
the second-order correction with additional fourth-order 
terms. These fourth-order corrections are generally 
smaller than the second-order ones, so they will only play 
a role when the second-order contributions are zero. 

We also find that E 2 is actually not always negative, 
and its sign will depend on the value of two competing 
terms, one of which is expressed in terms of (15) and the 
other one being 

S3 =yd"fe j Ki(fe, q)+6(fc, g)+6(fe, q)+U{k, g)], 

(73) 


where each of the ^i{K) are given in Appendix [C| 

It can be checked that, as a general trend, E 2 does 
not become negative when Ei is positive (i.e, when 
< 0 ), therefore the fourth-order analysis of entan¬ 
glement harvesting does not add significatively new re¬ 
sults. For illustration, we show in Fig. the two com¬ 
peting terms S 3 and {jCab^ in (72) for the case of 3+1 
dimensions with sudden and Gaussian switching. Recall 
that there would only be fourth-order contribution to en¬ 
tanglement if \Cab\‘^ > ^ 3 - 
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a) 


b) 


c) 





FIG. 2 . Dependence of entanglement harvesting on the size of the detectors for various values of the internal energy gap. The 
detectors are located in a 3+1-dimensional flat spacetime, and are switched on using the Gaussian prohle ( 22 ). We see that 
the larger the detectors, the less efficient they are for entanglement harvesting. 




FIG. 3. S 3 and \Cab\^ for the cases of spatially smeared 
detectors in 3+1 dimensions and a) Gaussian and b) sudden 
switching functions, for a time delay of A = lOT between their 
switchings. The vertical dashed lines represent the boundaries 
of the lightcone. Recall that there is a fourth-order contribu¬ 
tion to entanglement harvesting only if \Cab\^ > S 3 , which 
never happens in these cases. Notice the use of different y-axis 
scales for \Cab\^ and S 3 . 


IV. HARVESTING OF MUTUAL 
INFORMATION 

In the previous section we have shown that entangle¬ 
ment harvested between two two-level particle detectors 
strongly depends on the detectors’ switching, their spa¬ 
tial smearing, their internal energy gap and the dimen¬ 


sionality of spacetime. 

It is also an interesting question to find how particle 
detectors can harvest general correlations from the field, 
and not only entanglement. For instance, is it possible 
to harvest classical correlations from the field while the 
detectors remain spacelike separated?. 

In this section we will show that the harvesting of cor¬ 
relations from the field is much easier than harvesting 
entanglement. Namely, a pair of detectors prepared in 
their ground state can harvest correlations for a much 
wider range of scenarios than those that allow for entan¬ 
glement harvesting. 

We will characterize the total amount of correlations 
with the mutual information. 


Hpab) = S{pa) + S{pb) - S{pab), 


(74) 


where pA = Tr^(py^^) (and vice versa) is the partial 
trace of pab with respect to subsystem B {A) and S{p) = 
—Tr(/)logp) is the von Neumann entropy. 

For the state given by (13), the partial subsystems are 
described by the following density matrices 


so, once expanded to leading order in powers of the cou¬ 
pling strength, the mutual information takes the form 


Hpab) =C+ log(£+) + £_ log(£_) 

— CaA log(£AA) — ^BB log(£BB) + 

(76) 


where the quantities jC± are defined by 

>c± = 2 {^AA + ^BB ± \l{^AA ^BbY + 4 \Cab\^^ • 

(77) 

Note that in this expression the competition between 
the local and nonlocal (>C±) terms is also explicit. 

In the same fashion as in section m we show in Fig. 
1^ a collection of representative examples of the mutual 
information for all the cases defined in sections III Al and 

EH 
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FIG. 4. Mutual Information for all the cases described in sections [lI A| and |ll B[ in the following order: (a) 3D Gaussian switching 
and Gaussian spatial prohle; (b) 3D Gaussian switching and point like detectors; (c) 3D sudden switching and Gaussian spatial 
profile; (d) 3D sudden switching and near-pointlike detectors (cr/T 1); (e) ID Gaussian switching and Gaussian spatial 
profile; (f) ID Gaussian switching and pointlike detectors; (g) ID sudden switching and Gaussian spatial prohle; (h) ID sudden 
switching and point like detectors. The graphs show the harvested mutual information as a function of the spatial separation 
of the detectors and the time delay between their switching functions for flT = 1. In all plots the dashed lines represent the 
boundaries of the lightcone. These boundaries are placed at A = d ± T in the cases of point like detectors which are switched 
suddenly, and at A = d± 7T/v^ in the rest, which is a reasonable estimation given the discussion in section II, subsection Al. 
All cases of Gaussian spatial prohle have cr/T = 1 except for the near-pointlike detectors (plot d), for which a jT — 0.01. 


As we expected from the discussion in the study of 
the negativity, the mutual information maximizes inside 
the lightcone. In the same way as for the negativity, the 
total correlation between the detectors peaks inside the 
region of light contact, where direct exchange of real held 
quanta is possible. 

We hrst note that, unlike entanglement, harvesting of 
mutual information is possible for the whole spacetime, 
even for large temporal and spatial separations between 
the detectors long after entanglement harvesting is no 
longer possible (see Fig. |^. 

Due to this, when the separation between the detec¬ 
tors is large enough, we can identify the main source of 
harvested correlations as classical correlations present in 
the held vacuum. Nevertheless, we should not forget that 
other sources of quantum correlations (quantum discord 
[32]) may be present, and their study in future works may 
be of interest. 

In Fig. [^ we see how the dimensionality of space- 
time inhuences the harvesting of correlations: the gain 
of mutual information during light contact is more efh- 
cient in 3+1 dimensions than in 1+1. Interestingly, for a 
given value of QT the amount of correlations harvested 
in the spacelike-separation region is higher in the 1+1- 
dimensional case. 

We also observe the same effect of “damping and leak¬ 
age” that was observed in section m when increasing 
the internal energy gap is observed in the case of corre¬ 
lations: increasing QT decreases the overall value of the 


mutual information but also makes its support outside 
the lightcone increase. 


OT= 1 , A/T = 2.1 



FIG. 5. Mutual information (dashed red) and negativity 
(solid blue) to 0{Xl) for 1 + 1 -dimensional detectors with 
Gaussian spatial smearing and Gaussian switching. The 
black, dashed line represents the usual estimation of the ef¬ 
fective end of the lightcone /I = 7 + 7/\/2. 


V. CONCLUSIONS 

We have performed a detailed study of the phe¬ 
nomenon of entanglement and mutual information har- 
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vesting from the vacuum state of a scalar field, using a 
pair of Unruh-DeWitt particle detectors. 

First, we have analyzed the influence on this phe¬ 
nomenon of various aspects, such as the dimensionality 
of spacetime, the nature of the switching of the detectors 
(sudden vs. smooth), the detectors’ physical size and 
their internal energy structure. 

While the dimensionality of spacetime does not seem 
to strongly influence entanglement harvesting, the gain of 
mutual information during light contact is more efficient 
in 3-1-1 dimensions than in 1+1 dimensions (e.g, inside 
an optical fiber). Conversely, the amount of spacelike- 
correlation harvesting is noticeably larger in the 1+1- 
dimensional case. 

We have made a comparative study of two different 
kinds of switching: I) smooth Gaussian and II) sudden 
switching. We have shown that the smooth Gaussian 
switching is much more efficient than the sudden switch¬ 
ing to harvest entanglement from the vacuum in all cases: 
spacelike, timelike and lightlike. Remarkably, we have 
found that it is not possible to harvest spacelike entan¬ 
glement with sudden switching, neither in 3+1 nor 1+1 
dimensions, for the parameter landscape studied (with 
the exception of point like detectors in 1+1 dimensions for 
a very limited range of detector configurations and small 
spatial separations very close to the lightcone). This diffi¬ 
culty of harvesting entanglement with abruptly switched 
detectors contrasts with the case of smooth Gaussian 
switching, where it is always possible to find detector 
configurations for which there is significant entanglement 
harvesting for spacelike separated detectors arbitrarily 
far away from the lightcone. This result is striking, and 
runs perhaps contrary to previous conjectures (suggested 
by results obtained with non-perturbative harmonic os¬ 
cillator detectors [26]) that sudden switching may aid 
in generating entanglement. Our result is especially in¬ 
teresting in the light of the pioneering studies on en¬ 
tanglement harvesting that showed that super-oscillatory 
switching functions can enhance the detectors’ ability to 
harvest entanglement [T2|. 

However, the fact that Gaussian switching is more ef¬ 
ficient than sudden switching to harvest entanglement 
stems from the following reason: the more sudden the 
switching is, the stronger the local noise becomes 

[28] . The local noise competes with the nonlocal terms 
(Al) which give rise to entanglement. The smooth Gaus¬ 
sian switching strongly attenuates this noise and there¬ 
fore facilitates that the nonlocal terms dominate over the 
local noise terms. 

As for the size of the detectors, we have shown that 
if the detectors are smaller than the characteristic inter¬ 
action timescale, their harvesting abilities do not sub¬ 
stantially deviate from the pointlike case. When the size 
of the detectors becomes comparable to the interaction 
timescale, their ability to harvest entanglement gets in¬ 
creasingly hindered as their size increases. 

The energy gap of the detectors has a strong influ¬ 
ence in the extraction of entanglement and correlations 


from the vacuum, but this influence is radically different 
for smoothly switched detectors and suddenly switched 
detectors. For Gaussian switching it is possible to find 
values of the gap for which there is always spacelike en¬ 
tanglement harvesting. In more detail, increasing the di¬ 
mensionless parameter VtT leads to a “damping and leak¬ 
age” effect, which allows entanglement to be harvested in 
a broader range of spatiotemporal distances at the cost 
of harvesting less entanglement. This is not the case for 
sudden switching, where increasing the energy gap does 
not necessarily mean an increase in the area where en¬ 
tanglement can be harvested. This is consistent with our 
previous result that Gaussian switching is generally bet¬ 
ter to harvest vacuum entanglement. 

We have also discussed that harvesting classical corre¬ 
lations (and possibly quantum discord) is generally easier 
than harvesting entanglement. In particular it is pos¬ 
sible to harvest mutual information when we place the 
two detectors anywhere in the whole spacetime. Even 
though the amount of harvested correlation quickly de¬ 
creases with the spatiotemporal distance, it only vanishes 
in the limit of infinite separation. 

Finally, we have analyzed the first subleading-order 
correction to the negativity. We have seen that the next 
order contribution to the negativity vanishes generally at 
smaller spacetime distances than the second-order contri¬ 
bution, and thus going beyond leading order does not re¬ 
veal new phenomenology. Therefore a leading-order per¬ 
turbative approximation is generally enough to identify 
the regimes in which Unruh-DeWitt detectors can har¬ 
vest quantum entanglement from the field vacuum. 


ACKNOWLEDGMENTS 

The authors gratefully acknowledge Luis J. Garay for 
his valuable comments on the final version of this paper. 
We also thank Petar Simidzija for his inestimable help 
reviewing the latest versions of this manuscript. The 
work of E. M.-M. is supported by the National Sciences 
and Engineering Research Council of Canada through the 
Discovery program. 


Appendix A: Explicit computation of G 2 (/^), and 
M in the 3D, Gaussian switching case. 

Obtaining a closed-form expression of G 2 {i^) requires 
a more subtle process than Gi(K,r^) due to the nesting 
of the Ti, r 2 integrals. Recall that we start from 

/ OO PTI 

OO J —OO 

X ) . (Al) 

The innermost integral can be evaluated straightfor- 
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wardly, yielding 


.V^^2 


drije 


-(ri-7)^-i'ri(|K|-a)-5(|K|+a)^ 


1 + erf ( n - -i(|K| +a) 


+ e 




1 + erf 


^ri- 7 -ii(|K|+a)^ |. (A2) 


Notice that (A 2 ) has four summands. Two of them 


are combinations of exponential and Gaussian functions, 
whereas the other two are of the form 


/ oo 

dy e~^ -^^y-y (^y _ ^ 

-oo 


(A3) 


Partial differentiation under the integral sign of the 
previous expression with respect to the parameter a 
yields 


d 

da 


I{a,b)= r 

J —( 


dy 


- 2 ae-“'erf (y - ia) 


-i^e 

Vtt 




= -2a/(a,6)-iy2e-5('’-2“)" 


(A4) 


This is a linear, first-order, nonhomogeneous diferen- 
tial equation, which can be solved, for instance, via varia¬ 
tion of constants. The solution to the homogeneous equa¬ 
tion is simply 


lH{a,b) = Ce “ , 


(A5) 


and now, allowing the constant to be a function of a and 
6 , and inserting b) in Eq. (|A4|), we obtain 


A 

da 


C{a,b)e-^^ = -i\/2e-5(^-2“)' 


C{a,b) = —i-y/Tre * erfi 


2(2 T b 

2 V 2 


So the closed form on the integral /(a, b) is 


I{a^b) = —iy/ne ^ 4 erfi 


V2 


(A 6 ) 


(A7) 




G2{k) 


Therefore, using this in ( |A 2 [ ) we get 

+ i {^(\ k \ +a),|K| - a+ 217) 


(A 8 ) 


Now that we have closed expressions for Gi(k, r^) and 
we can compute and M. Recall that, from 
(pl and ra. 


=Xyd-^k ^ ' 3 . ,,, gi(K, r^)Gl{K, T^), (A9) 

(27rj Z\k\ 

_j , /\ rj* _^ ^ 0- ^ 

CAB=X^ld^k ^ 'hi. gi(/c,0)gt(«,7), (AlO) 


,=A“/, 


\M\ =A2 


I 


(27r)^2|fe| 
d3fegife-(®-4-a!B) 


g-ife 2 <T = 


(27r)^2|fe 


G2 {k) 


■ (All) 


Writing these integrals in terms of the dimensionless 
variable k, = kT^ using spherical coordinates and choos¬ 
ing the 2 (-axis in the direction of k we obtain: 

\ 2 poo 


47r2r2 jQ 
'87r(l + (52) 






(A 12 ) 


which gives the values for Gaa and jCbb^ 
12 /»! 
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being jCba = and, finally, 

|A4| = 


(A13) 
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Appendix B: Explicit computation of S 2 {h 


To find simple expressions for 5 ' 2 (^): we will separate 
the cases in which the switching functions do and do 
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not overlap. In the nonoverlapping case, let us assume 
without loss of generality that detector A is switched off 
before detector B is switched on. In this case, the first 


contribution of (42) is trivially zero, and in the second 


one the time integrals decouple. Substituting the dimen¬ 
sionless variables = t^/T we get 


Appendix C: The fourth-order term S 3 

The four components ^^(^, 77 ), i = 1,..., 4 of S 3 are 

[Fimm? 


^i{k,q) =Xa^b 


^off POO nti nOO 

S'2non(«) dn / fiT^ei»{Ti+T 2 )^-i\K.\iTi-T 2 ) X dti dt2 dt[ dt 2 XA{ti)XBit 2 )XAit'i)XB{t 2 ) 

I on I on J —00 J —00 J —00 J —00 

F -D F A 


4|fc||q| 

^00 




— a‘^) T ^ 


(Bl) 


where 

When the switching functions overlap, a more careful 
analysis is required. We assume without loss of generality 
that detector A is switched on simultaneously, or before 
detector B. For illustration let us further assume that 
detector B is switched off simultaneously or after detector 
A. The total interaction time interval [T^^,Tg^] can be 
then subdivided into different regions: 

• Regions of no overlap: U [T^^, Tg^]. 

• Region of overlap: 

Thus we will have three different contributions, 

S2{k) = + 

(B2) 

I n th e nonoverlapping regions we get a result analogous 
to (BI) —the first contribution is zero and the integrals 


gia(ri+r2)g-i|K|(ri-r2) 


decouple in the second-, so we obtain 

^off ^on 

S2^^-Fb-'\k)=T^[ "dn [ "dr2 

J^on 

7^2gi7(a-|/*|) 

a‘^ — K? 

X _gi[«(7+l) + l«^l(7-l)] _^gi7(a+l«^l) ^ 




ia(ri+r 2 )g-i|/«|(Ti-r 2 ) 


2^2 ia(7+l),_ 

:I___jgiQ^(l-7) _^i(7- 

— tx? V 


l)kl _gi(«-7kl) _^gi(«7-kl)^ ^ 


(B3) 


In the overlapping region both contributions of (|^ 
are equal and nonzero, so we can write 


=2T^ rdT2e' 


ia(ri+r 2 )g-i|/«|(ri-r 2 ) 


= 2 T^ 


gi[|K;|(7-l)+a(7+l)] _ g2i7a g2i7a _ ^2\a 


2a{\K\ + a) _ 


(B4) 


Adding all three contributions, we obtain (45). 


X e 
X 


—00 ^—00 «/—00 ^—00 




gi|^l(^i-^i)gikl(^ 2 -t 2 ) 


(Cl) 


^ 2(^5 9 ) —Xa^b 


[F(fc)] 2 [F(q )]2 


4|fe||q| 

00 pti pOO pt'^ 


/ OO pZ\ pOO pZ^ 

dti / dt2 / d7 / d4 XA{tl)XB{t2)XB{t'i)XA{t2. 

-OO J —OO J —OO J —OO 




J\k\(ti-t[) ^i\q\(t2-t2) ^-iik-q)-(xA-XB) 
J\k\(t2-ti) i\q\(t[-t2) i(k-\-q)-{xA-XB) 


gi|fc|(ti-t 2 )gikl(^ 2 -ti) 


(C 2 ) 


6(fe,9) =XaXb 


[F(fe)] 2 [F(g )]2 


4|fc||qr| 

/ OO pti pOO pt^ 

dti / dt2 / dt[ / dt2 XB{ti)XA{t2)XA{t[)xB{t2] 

-OO J —OO J —OO J —OO 


^\\k\{ti-t'^) ^\\q\{t2-t'2) ^i{k-q)-{xA-XB) 


^i\k\{t2-ti) i\q\{t'^-t'2) -i{k+q)-{xA-XB) 


+ e 


i|fc|(ti-t2)gi|q|(t2-ti) 


(C3) 


U{k,q) =X\X'^ 


[F{k)nF{q)f 


A\k\\q\ 

/ OO pti pOO pt^ 

dti / dt2 / dt'i / d4 XB(ii)XA(i2)XB(ii)Xvi(i7 

-OO J —OO J —OO J —OO 


X 


gi|fc|(t2-ti)gi|q|(ti-t2)g-i(^-9)-(®A-®B) 


gi|^l(^i-^2)gikl(^2-ti)gi(fc-q)-(®A-®B) 


+ e- 


i|fc|(ti-ti) i|q|(t2-t2) 


(C4) 


To evaluate these integrals, let us consider the simpler 
case where the detectors are identical and the detectors’ 
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switching functions have supports that are far apart. In 
this case the only nonzero contribution to S 3 is ^ 4 , which 
can, under these assumptions, be split into three parts 
that we denote as: 

Uaik,q) 

^ (C5) 


Ubik, q) r,), 

^ (C 6 ) 


Uc{k,q) =A 


,[F{k)nF{q)Y 


4\k\\q\ 


R3{k,v), (C7) 


where k = kT and rj = qT are dimensionless momenta 
and the functions q) are given by 


/ OO nOO nOO nOO 

dri / dT 2 / dT{ / drg 

00 J —00 J —00 J —00 


— 00 J—00 J—00 J—00 

Xa{ti) Xb {t2)xa{t[)xb (t2 ) 

Ja{ri-\-r 2 -T[-T^) i\K\{r 2 -ri) i\r]\{r[-T^) 


/ OO pOO pOO pOO 

dri / dT 2 / d 4 / dT 2 

-OO J —OO J —OO J —OO 

Xa{ti)xb{t2)xa{t[)xb{t2) 

^ gia(ri+r2-Tl-T2)gi|/«|(ri-T2)gi|T7|(r2-rl) 

(C9) 

/ OO pOO pOO pOO 

dri / dr 2 / dr^ / dr^ 

-OO J —OO J —OO J —OO 

XtI (n )XB (r2) XA (rO XB (r2) 

^ gia(ri+T2-ri-r2)gi|/«|(ri-Tl)gi|T7|(r2-r2)^ 

(CIO) 

Thus S 3 can be written as 

poo poo 

E3={4nfT^X^ d\K\ d\rj\[F{k)nF{q)Y 

Jo Jo 

sin(/3|K|) sin(;d|r 7 |) _ , , ^ . 

77 ) + i?2(«, 77 )] 


(Cll) 


|K||r7|i?3(«,^?)|- 


(C 8 ) 


In the case when the detectors’ switching functions are 
two Gaussians far apart enough so that their mutual 
overlap can be neglected, all the previous integrals ad¬ 
mit closed-form expressions. Using (20) and (22) we get 
to the final expression 


-3 =- 


X^e 1 + 5 ^ 
1287r/32 (1 + (^ 2 ) 

(/3-7)^ 


_ (P + o)^ 

g 2(1 + 52) 


1 +ierfi 


( /3 + 7 \ 

\V 2 VTTP). 




g 2(1+52) 


1 —isign(/3— 7 )erfi 






1 — i erfi 


(/3 + 7)^ 


_g2(l + 52) 


1 -h i sign (/d — 7 ) erfi 


1287r/32 (1 + 52 ) 

[a-i(^ + 7 )]^ 

e 2(1+52) 


P + X 

erfc - i(/? - 7)^ (a+ i{pFx) 


l/?- 7 l 

\/ 2 \/rT^ 


)]} 


\p-i\ 


V 2 V 1 + 52 




Q-i(^ + 7)]^ / Q/ — + 7 ) ^ [Q + i(/3-7)] 


V \/2\/rTF 


VV 2 ^/^+^y ® 

-e‘7w+erfc7 + i(/3-7) 


V ^/2^/^T^ 


A+ 


647r2 (1 + 52 ) 


2\/1 + 52 — \/^Q:e2(i+'5^) erfc 


In the case of nonoverlapping sudden switchings, the different Ri{K,r]) read 

p-i[2a-(7-l)|r7| + (27+l)|K,|] 

Ri{fi,v) = 

-R 2 («,■??) = 


(C 12 ) 


(^_gi« _|_ gi(2a+lnl) _ gi(«-l-2|nl) _|_ g'lul^ ^gi(«-l-kl) _ (^gi(a-l-7|K'l) _ gi(7-l-l)|K|^ 


R3{h,v) =- 


(^Q ;2 _ (q ;2 _ 1 ^^ ^ 

g-i[2a+(7+l)|r7|-(7-l)|K|] \2 / \ 

X_ ( gi(«-l-|»7l) _ 1 \ (gi(a+|K|) _ 1 \ 

(a+|77|)2(a+|K|)2 I M ; 

sin 2 sin 2 (^ 7 ^) 

(a+|r7|)2(a+|«|)2 ’ 
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(C13) 
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and 


Z 3 will be given by (Cll), upon substitution of the Ri{K^r]) integrals. 
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